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ABSTRACT 

The  Euler  decomposition,  when  applied  to  the  polarization  scattering  matrix,  attempts  to  extract 
phenomenological  information  about  the  scattering  target.  Because  the  Euler  parameters  constitute  a  more  physically 
relevant  set  of  parameters  than  the  traditional  HE1-VV  ISAR  representations,  they  have  potential  to  improve  ATR 
performance.  The  Euler  parameter's  usefulness  in  target  recognition,  however,  is  effected  by  several  layers  of  signature 
variability.  Unfortunately,  many  of  the  variability  layers  are  often  omitted  in  a  typical  ATR  study.  A  complete  ATR 
algorithm  was  therefore  developed  that  allows  for  all  layers  of  variability  and  requires  no  previous  knowledge  of  the 
target's  position,  orientation,  or  average  reflectivity.  The  complete  ATR  algorithm  was  then  used  to  assess  the 
effectiveness  of  Euler  ISAR  imagery  in  target  recognition  when  all  layers  of  variability  are  considered.  The  general 
approach  and  sub-methods  used  to  construct  the  complete  ATR  system  will  be  presented,  including  the  methods  to 
determine  the  targets  orientation,  registration,  and  to  compare  it  to  a  library  of  pre-rendered  target  images.  Finally,  the 
performance  of  the  Euler  parameters  in  target  recognition  using  the  complete  ATR  algorithm  will  be  presented. 
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1.  INTRODUCTION 

When  a  polarimetric  radar  system  is  used  to  measure  the  electromagnetic  backscattering  of  an  object  at  each 
possible  polarization,  the  complex  Radar  Cross  Section  (RCS)  values  can  be  formed  into  the  polarization  scattering 
matrix  known  as  the  Sinclair  matrix  [see  Eq.  (1)]. 
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In  this  way,  the  Sinclair  matrix  S  links  the  transmitted  fields  Ex  to  the  received  electric  fields  Er  for  the  traditional 
polarizations  of  Horizontal  (H)  and  Vertical  (V).  The  wave's  propagation  term  in  parentheses  in  Eq.  (1)  is  factored  and 
normalized  out. 


Use  of  a  polarimetric  Inverse  Synthetic  Aperture  Radar  (ISAR)  system  enables  the  formation  of  an  object's 
spatial  two-dimensional  image.  The  ISAR  image  contains  a  spatial  map  of  the  object's  scattering  properties  so  that  the 
Sinclair  matrix  is  known  for  each  pixel  in  the  image.  This  image  formation  is  possible  because  scattering  measurements 
made  for  a  sweep  in  frequencies  as  well  as  for  a  sweep  in  azimuthal  look  angles  are  spectral  assessments  of  the  object's 
spatial  extents.  When  such  a  sweep  of  RCS  values  for  frequency  and  azimuth  are  Fourier  transformed,  the  spatial 
distribution  of  the  object's  scattering  properties  are  recovered  in  a  pixelized  image. 
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The  Sinclair  matrix  contains  all  of  the  scattering  information  available  about  the  objects  contained  in  the 
corresponding  pixel,  but  is  in  a  form  that  is  not  particularly  intuitive  or  phenomenological.  There  exists  various 
decomposition  techniques  that  attempt  to  transform  the  Sinclair  matrix  into  parameters  that  have  more 
phenomenological  meaning1.  In  particular,  the  Euler  decomposition  is  one  that  seeks  to  diagonalize  the  Sinclair  matrix 
through  a  coordinate  rotation  and  then  extract  the  Euler  parameters  as  defined  by  Kennaugh  and  Huynen2'3.  When 
applied  to  the  scattering  measurements,  the  Euler  decomposition  leads  to  ISAR  images  of  the  object's  Euler  scattering 
properties. 

A  conjugate-similarity  unitary  transform  [see  Eq.  (2)]  is  the  operation  required  to  successfully  diagonalize  the 
Sinclair  matrix  S,  as  demonstrated  first  by  Autonne4.  It  should  be  noted  that  the  Sinclair  matrix  is  defined  in  the  Back 
Scattering  Alignment  (BSA)  coordinate  system,  in  contrast  to  the  Forward  Scattering  Alignment  (FSA)  used  in  optics. 

Sd=UtSU  (2) 


As  a  result  of  the  BSA,  the  conjugate  eigenvalue  equation  must  be  solved  to  construct  the  transform  matrix  U,  instead  of 
the  traditional  eigenvalue  equation  used  in  optics.  The  Euler  parameters  are  then  defined  in  terms  of  the  transform 
matrix  U  as  shown  in  Eq.  (3)  and  the  diagonalized  scattering  matrix  Sd  as  shown  in  Eq.  (4). 
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These  five  Euler  parameters  describe  the  scattering  object's  maximum  reflectivity  m,  orientation  angle  i//.  symmetry  r, 
bounce  angle  v,  and  polarizability  y.  The  physical  meaning  of  the  Euler  parameters  can  be  understood  by  illustrating  the 
parameters'  possible  values  and  some  simple  shapes  representative  of  these  values,  as  shown  in  Fig.  1. 

To  use  the  Euler  parameters,  the  definitions  in  Eqs.  (3)  and  (4)  must  be  inverted  to  yield  transform  equations 
establishing  the  Euler  parameters  as  functions  of  the  measured  Sinclair  matrix  elements.  This  inversion  is  best 
accomplished  by  switching  to  a  power  matrix  form  known  as  the  Kennaugh  matrix  that  is  equivalent  to  the  electric-field 
scattering  matrix.  The  Kennaugh  matrix  can  than  be  diagonalized  one  matrix  element  at  a  time  through  similarity 
rotations  similar  to  Eq.  (2),  and  the  Euler  parameters  are  extracted  one  by  one.  The  detailed  derivation  of  the  Euler 
transform  equations  has  been  presented  previously5,6. 

Once  the  Euler  transform  equations  were  derived  and  implemented,  the  test  object  Slicy  was  imaged  in  a 
compact  radar  range  as  an  initial  check  and  intuitive  demonstration  of  the  Euler  parameters.  As  found  in  Fig.  2,  Slicy 
contains  a  trihedral  [see  Fig.  2(a)]  that  shows  up  properly  as  having  an  odd  bounce  and  a  dihedral  [see  Fig.  2(b)]  that 
correctly  shows  up  as  even  bounce.  Slicy  also  has  back  corners  [see  Fig.  2(f)]  and  a  front  edge  [see  Fig.  2(d)]  that 
correctly  appear  as  polarizing.  It  is  interesting  to  note  that  Slicy  appears  almost  entirely  symmetrical  at  this  resolution. 
Upon  inspection,  Slicy  can  be  seen  to  contain  several  other  notable  scattering  features  that  show  up  in  its  Euler  ISAR 
images. 


Due  to  their  phenomenological  nature,  the  Euler  parameters  have  the  potential  to  improve  Automatic  Target 
Recognition  (ATR).  A  compete  ATR  algorithm  was  developed  that  includes  all  of  the  layers  of  target  variability 
encountered  in  the  field  and  requires  no  human  intervention  or  previous  knowledge  of  the  unknown  target.  The  complete 
ATR  algorithm  uses  the  Euler  ISAR  images  of  an  unknown  object  to  attempt  to  match  it  to  an  object  in  the  reference 
library.  A  persistence  optimization  method  was  also  developed  as  part  of  the  complete  ATR  algorithm  that  was  found  to 
improve  performance.  The  ATR  algorithm  was  tested  using  ISAR  signatures  of  tanks  obtained  in  previous  studies  under 
the  support  and  direction  of  the  U.S.  Army's  National  Ground  Intelligence  Center  (NGIC)7  8. 
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Fig.  1.  Sample  scattering  objects  that  display  the  extremes  of  each  Euler  parameter. 
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Fig.  2.  The  test  object  Slicy  and  its  associated  Euler  ISAR  images  measured  at  0°  azimuth  and  5°  elevation,  demonstrating 
the  physical  meaning  of  each  Euler  parameter  and  their  applicability  to  actual  measured  data. 


2.  IMPLEMENTATION  OF  THE  ATR  ALGORITHM 

The  ATR  algorithm  uses  enough  scattering  data  to  form  one  polarimetric  ISAR  image  of  the  unknown  object, 
attempts  to  match  its  image  to  a  reference  library  of  previously  acquired  images,  and  then  assigns  the  identification  of  the 
best  matching  library  object.  In  order  to  determine  the  most  effective  ATR  configuration,  the  algorithm  was  designed  to 


construct  and  use  images  in  any  of  the  traditional  parameter  spaces  HH,  HV,  VH,  VV,  as  well  as  in  any  of  the  Euler 
parameter  spaces  m,  y,  i//,  r,  or  v.  A  persistence  optimization  was  also  developed  that  improves  the  ATR  performance  of 
the  Euler  parameters. 

As  outlined  in  Fig.  3,  several  steps  are  involved  in  the  complete  ATR  algorithm.  The  library  of  reference  images 
must  be  formed  before  the  algorithm  can  be  used.  Each  object  to  be  included  in  the  library  is  imaged  in  a  radar  range  at 
every  possible  azimuthal  look  angle  to  within  1°  and  at  some  fixed  elevation  angle.  Each  library  image  is  transformed 
into  Euler  ISAR  images,  thresholded  to  reduce  noise,  exactly  back-rotated  to  the  common  axis  of  0°  azimuth,  and 
analyzed  to  determine  the  persistence-optimization  values.  The  azimuthal  persistence  of  each  pixel  is  calculated  as  the 
extent  in  the  azimuthal  sequence  of  ISAR  images  for  which  the  pixel's  values  do  not  change  beyond  some  tolerance 
band.  The  more  persistent  pixels  are  then  assigned  greater  weights  being  deemed  the  more  reliable  pixels.  The 
persistence  weights  are  stored  as  part  of  the  reference  library  and  are  used  during  image  comparisons  to  optimize  ATR 
performance. 

For  the  particular  library  suite  of  tanks  used  in  this  study,  a  522  GHz  compact  range  was  used  with  1/1 6th  scale 
replicas  imaged  at  a  5.00°  elevation  angle  to  match  a  corresponding  32.625  GHz  radar  in  the  field.  The  scaled  tanks  were 
fabricated  with  scaled  material  properties  through  the  ERADS  partnership  and  imaged  at  the  University  of 
Massachusetts  Lowell  Submillimeter  Technology  Laboratory  (STL)  under  the  support  of  the  U.S.  Army's  National 
Ground  Intelligence  Center  (NGIC). 

When  the  ATR  algorithm  is  ready  to  be  used,  the  scattering  data  of  the  unknown  data  is  first  analyzed  to 
establish  its  approximate  center  and  azimuthal  orientation  (see  Fig.  3).  The  object's  geometric  center  is  calculated  as  the 
center  of  the  distribution  of  reflected  power. 

The  approximate  azimuth  of  the  unknown  vehicle  is  found  by  assuming  a  rectangular  shape  and  then  finding  the 
angle  at  which  the  rotation  of  the  image  causes  the  vehicle's  rectangular  outline  to  become  parallel  to  the  image's  axes. 
This  occurs  at  the  angle  at  which  the  one-dimensional  profile  of  the  object  becomes  narrowest,  as  projected  on  the 
image's  x-axis.  Because  the  azimuth  cannot  be  determined  to  a  greater  accuracy  than  a  few  degrees,  several  possible 
azimuths  are  reported  and  imaged  for  the  one  unknown  object.  When  the  comparisons  are  later  made  for  each  possible 
azimuth,  the  azimuth  for  which  the  unknown  object's  image  best  matches  the  library  images  is  retained  as  the  true 
azimuth. 


Once  the  approximate  center  and  possible  azimuths  are  determined,  the  scattering  data  is  formed  into  the  ISAR 
images  using  Fourier  transforms.  The  image  is  back-rotated  as  part  of  the  image  formation  process  by  the  amount  of  its 
azimuth  angle  in  order  to  align  all  of  the  possible  images  of  unknown  targets  to  the  common  axis  at  0°  azimuth  for  proper 
comparison  to  the  library  images  which  where  already  back-rotated.  Instead  of  a  traditional  image  rotation,  a  more  exact 
back-rotation  was  developed.  The  traditional  image  rotation  suffers  degradation  because  pixels  are  assigned  to  the  new 
non-integer  rotated  locations  but  must  be  rounded  to  the  nearest  integer  pixel  location.  The  new  exact  rotation  method 
applies  the  coordinate  rotation  before  the  image  is  formed  by  the  Fourier  transforms,  and  thus  avoids  the  degradation  due 
to  rounding.  In  practice,  the  novel  rotation  technique  calculates  the  new  non-integer  rotated  locations  and  then  shifts  the 
Fourier  transform  to  those  locations  when  it  is  applied  as  shown  in  Eq.  (5). 
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Next,  the  Euler  transform  equations  are  applied  to  each  back-rotated  ISAR  image  of  the  unknown  object  in 
order  to  create  Euler  ISAR  images.  The  images  are  thresholded  to  reduce  noise  and  are  ready  to  be  compared  to  the 
library  images 

Each  of  the  unknown  object's  images  formed  at  the  possible  azimuths  is  compared  in  a  brute-force  fashion  to 
every  possible  library  object  at  that  azimuth.  The  comparison  of  two  ISAR  images  involves  registering  them  through  an 
autocorrelation  technique  and  then  computing  the  Average  Percent  Difference  (APD)  between  the  images.  The  APD 
method  of  establishing  a  correlation  score  has  been  presented  and  analyzed  elsewhere9.  While  the  individual  pixel 


percent  differences  are  being  averaged  to  form  the  APD,  the  persistence  weights  calculated  from  the  azimuthal 
persistence  of  the  library  images  are  applied  to  optimize  performance.  The  APD  scores  of  all  of  the  comparisons  are 
tabulated  by  the  algorithm  and  the  comparison  with  the  lowest  score  is  assigned  as  the  unknown  object's  identification. 


Fig.  3.  Schematic  process  flow  of  the  full  ATR  algorithm. 


3.  TESTING  THE  ATR  ALGORITHM 


A  separate  test  module  was  developed  which  runs  hundreds  of  unknown  images  from  a  test  suite  through  the 
ATR  algorithm,  compares  the  identifications  to  ground  truth  and  establishes  performance  curves.  All  of  the  APD 
comparisons  scores  are  grouped  into  histograms  and  normalized  to  establish  probability  density  curves.  The  probability 
density  curves  are  then  scanned  as  shown  in  Fig.  4  to  construct  Receiver  Operating  Characteristic  (ROC)  curves.  ROC 
curves  are  formed  by  plotting  the  area  of  true-positive  identifications  versus  the  area  of  false-positive  identification  for 
each  possible  decision  threshold.  The  closer  the  ROC  curve  approaches  the  top-left  corner  of  the  graph,  the  better  the 
ATR  performance.  Each  ROC  curve  therefore  represents  the  performance  for  a  specific  ATR  algorithm  configuration 
when  hundreds  of  tests  are  done  from  the  test  suite. 


The  performance  tests  were  carried  out  on  a  suite  of  similar  tanks  with  differing  equipment  configurations.  It 
must  be  stressed  that  the  task  of  differentiating  similar  tanks  that  differ  mainly  in  equipment  and  not  in  overall  size  and 
shape  is  among  the  most  difficult  tasks  that  an  ATR  algorithm  can  encounter.  This  challenging  test  suite  was  chosen  so 
as  to  best  reveal  the  effect  of  the  Euler  parameters.  Two  of  the  tanks  from  the  test  suite  are  presented  in  Figs.  5  and  6  to 
illustrate  how  similar  the  objects  are  that  must  be  differentiated  by  the  ATR  algorithm.  The  suite  of  library  images  was 


constructed  using  different  measurements  of  different  physical  copies  of  the  same  tank  configurations  as  in  the  test  suite. 
Each  of  the  tanks  in  the  test  suite  and  the  library  suite  was  measured  as  1/1 6th  scaled  model  replicas  in  the  522  GHz 
compact  radar  range  at  the  University  of  Massachusetts  Lowell  Submillimeter  Technology  Laboratory  (STL)  to  match 
an  equivalent  full-scale  radar  system  at  32.625  GHz.  The  scaled  tanks  were  built  through  the  ERADS  program  under 
sponsorship  of  the  United  States  Army's  National  Ground  Intelligence  Center  (NGIC)  as  part  of  previous  studies7  8.  Each 
tank  was  imaged  in  free  space  at  a  5.00°  elevation  angle. 


4.  ATR  PERFORMANCE  RESULTS 

The  ATR  algorithm  was  tested  using  the  entire  test  suite  for  different  imaging  parameters.  The  results  presented 
in  Fig.  7  indicate  that  the  angular  Euler  parameters  ///,  r,  v,  and  y  give  ROC  curves  further  away  from  the  top-left  corner 
and  thus  lead  to  worse  target  recognition  than  the  traditional  imaging  parameters  HH,  HV,  VH,  and  VV.  This  indicates 
that  although  the  angular  Euler  parameters  are  more  phenomenological  than  the  traditional  parameters,  they  are  sensitive 
enough  to  noise  to  yield  worse  performance. 

Adjusting  the  scales  on  Fig.  7  to  zoom  in  on  the  top-left  corner  reveals  in  Fig.  8  that  the  Euler  parameter  m  leads 
to  the  best  ATR  performance,  better  than  the  traditional  parameters. 
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Fig.  7.  ROC  curves  showing  ATR  performance  for  the  various 
imaging  parameters  without  persistence  optimization. 
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Fig.  8.  Zoomed-in  ROC  curves  of  ATR  performance  for  the  various 
imaging  parameters  without  persistence  optimization. 


The  results  in  Figs.  7  and  8  were  obtained  when  no  azimuthal  persistence  optimization  was  used.  When 
persistence  optimization  was  turned  on,  the  ATR  algorithm's  performance  improved  for  the  angular  Euler  parameter 
configurations  but  did  not  change  much  for  the  magnitude  parameter  configurations,  which  includes  the  Euler  parameter 
m  and  the  traditional  parameters  TUT,  E1V,  VH,  VV.  Only  the  ATR  performance  results  for  the  persistence-optimized 
Euler  polarizability  parameter  and  the  traditional  parameter  HV  are  shown  in  Figs.  9  and  10,  but  they  represent  all  of  the 
angular  parameters  and  magnitude  parameter  respectively. 

The  reader  must  remember  that  each  ROC  curve  represents  the  total  target  recognition  performance  for  a 
specific  ATR  algorithm  configuration  when  all  of  the  tanks  in  the  suite  are  tested  at  every  possible  angle,  and  thus  each 
curve  has  broader  significance  then  any  one  individual  comparison. 

Even  with  the  improvements  gained  using  persistence  optimization,  the  angular  Euler  parameters  still  yield 
worse  target  recognition  than  the  traditional  parameters.  It  would  seem  that  aside  from  the  Euler  parameter  m,  the  Euler 
parameters  are  useless  in  improving  target  recognition.  However,  the  parameters  have  only  been  treated  separately  up  to 
this  point.  Because  the  Euler  parameters  represent  independent  scattering  properties  of  the  exact  same  object,  they  can 
still  improve  target  recognition  when  used  in  a  consolidated  manner.  When  implementing  consolidation,  the  scattering 
data  is  still  imaged  into  independent  Euler  images  and  the  different  Euler  parameter  images  are  still  compared  to  the 
library  images  independently,  but  the  independent  correlation  scores  are  consolidated  to  improve  performance. 


probability  of  false-positives  identification  (%) 

Fig.  9.  ROC  curves  of  ATR  performance  using  polarizability 
images  with  and  without  persistence-optimization, 
representative  of  all  of  the  angular  Euler  parameters. 
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Fig.  10.  ROC  curves  of  ATR  performance  using  HV  images 
with  and  without  persistence-optimization,  representative  of 
all  of  the  magnitude  parameters  HH,  HV,  VH,  VV,  and  m. 


The  persistence-optimized  Euler  parameters  were  combined  into  a  consolidated  score  s  according  to  Eq.  (6).  In 
this  equation,  sm,  for  example,  is  the  correlation  score  between  the  unknown  vehicle  and  the  library  vehicle  when  they  are 
both  imaged  in  Euler  parameter  m  space.  The  centering  terms  ,v,„o,  syo,  srf),  svth  s„o,  in  Eq.  (6)  are  fixed  constants  that  are 
acquired  by  optimizing  them  for  peak  ATR  performance  using  a  training  set  of  vehicles  separate  from  the  test  suite. 


Conceptually,  the  five  independent  Euler  parameter  correlation  scores  specify  a  point  in  five-dimensional  space. 
The  consolidated  score  s  defined  in  Eq.  6  is  then  the  distance  in  this  five-dimensional  space  from  the  correlation  point  to 
a  fixed  reference  point.  The  shorter  the  distance  s,  the  better  the  two  images  match. 

The  final  ATR  performance  results  are  presented  as  zoomed-in  ROC  curves  in  Fig.  11.  The  Euler  parameter  m 
yields  the  best  target  recognition  when  treated  separately,  even  better  than  the  traditional  HH  imaging  space.  However, 
the  best  ATR  performance  was  achieved  by  consolidating  the  persistence-optimized  Euler  parameters  into  a  single 
correlation  score  s.  It  should  be  remembered  that  these  ROC  curves  represent  the  target  recognition  performance  of  the 
ATR  algorithm  for  the  challenging  test  suite  of  spatially  similar  tanks  differing  mostly  on  the  level  of  equipment. 
Inclusion  of  spatially  different  tanks  such  as  the  Ml  in  the  test  suite  were  found  to  only  improve  the  ATR  performance 
results  and  push  the  ROC  curves  further  into  the  top-left  corner. 

The  ROC  curve  midpoint,  shown  as  a  box  in  Fig.  11,  indicates  that  the  ATR  algorithm  using  HH  images  leads 
to  a  98.9%  true-positive  identification  rate.  In  contrast,  the  Euler  parameter  m  gives  an  ATR  performance  curve  midpoint 
at  99.4%  true-positives.  Best  of  all,  the  ATR  algorithm  based  on  the  consolidated  Euler  score  s'  gives  a  ROC  curve 
midpoint  with  a  99.7%  true-positive  identification  rate. 


Fig.  11.  Zoomed-in  ROC  curves  showing  ATR  performance  for  the  best  three  imaging  configurations. 


5.  CONCLUSIONS 

A  new,  full  ATR  algorithm  has  been  developed  using  a  reference  library  of  high-fidelity  radar  images,  a  novel 
azimuth-identification  method,  and  novel  imaging  techniques  which  is  shown  to  successfully  recognize  targets  at  high 
rates.  Setting  the  ATR  algorithm  to  use  persistence-optimized  Euler  parameter  images  with  consolidated  correlation 
scores  resulted  in  better  target  recognition  than  the  traditional  parameter  approaches. 
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